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Using a (2+l)-d viscous hydrodynamical model, we study tlie dependence of fiow observables on 
the collision energy ranging from -^s = 7.7AGeV at the Relativistic Heavy Ion Collider (RHIC) 
to -^s = 2760 A GeV at the Large Hadron Collider (LHC). With a realistic equation of state, 
Glauber model initial conditions and a small specific shear viscosity rj/s = 0.08, the differential 
charged hadron elliptic fiow V2^{pt,^) is found to exhibit a very broad maximum as a function 
of ^/s around top RHIC energy, rendering it almost independent of collision energy for 39 < \/s < 
2760 A GeV. Compared to ideal fluid dynamical simulations, this "saturation" of elliptic fiow is 
shifted to higher collision energies by shear viscous effects. For color-glass motivated MC-KLN 
initial conditions, which require a larger shear viscosity rj/s = 0.2 to reproduce the measured 
elliptic flow, a similar "saturation" is not observed up to LHC energies, except for very low pr- We 
emphasize that this "saturation" of the elliptic flow is not associated with the QCD phase transition, 
but arises from the interplay between radial and elliptic flow which shifts with depending on the 
fluid's viscosity and leads to a subtle cancellation between increasing contributions from light and 
decreasing contributions from heavy particles to V2 in the ^/s range where vf^{pT,^) at flxed pr 
is maximal. By generalizing the definition of spatial eccentricity to isothermal hyper-surfaces, we 
calculate on the kinetic freeze-out surface at different collision energies. Up to top RHIC energy, 
^/s = 200y4GeV, the fireball is still out-of-plane deformed at freeze out, while at LHC energy the 
final spatial eccentricity is predicted to approach zero. 



I. INTRODUCTION 

The first decade of experiments at the Relativistic 
Heavy Ion Collider (RHIC) has established the existence 
of a strongly coupled Quark-Gluon Plasma (sQGP), a 
new state of nuclear matter with partonic degrees of free- 
dom. One of the major goals of heavy-ion collision exper- 
iments right now is to explore the QCD phase diagram. 
The recent Beam Energy Scan (BES) program [iHal at 
RHIC is motivated by searching for the phase bound- 
ary between normal nuclear matter and sQGP as well as 
for the theoretically predicted QCD critical point [7l-[l0|. 
The BES program at RHIC together with Pb-hPb cob 
lisions at the Large Hadron Collider (LHC) provide us 
with a unique opportunity to study systematically the 
collision energy dependence of relativistic heavy-ion col- 
lision observables. 

The study of collective flow in relativistic heavy ion 
collisions has the potential to offer insights into proper- 
ties of the produced matter. Anisotropic flows, especially 
the elliptic flow V2, are widely studied in heavy- ion colli- 
sions. In the mid-rapidity region, the dependence of ellip- 
tic flow on transverse momentum and collision energy are 
crucial for our understanding of the properties of sQGP. 
Elliptic flow can provide information about the specific 
shear viscosity, which controls the conversion efficiency 
of anisotropic spatial pressure gradients to momentum 
anisotropics in a hydrodynamic description, and about 
the equation of state of the matter created at early times 
[lll428l | . Recent measurements of higher order anisotropic 
fiow coefficients, w„(n > 3) have generated strong inter- 
est due to their abili ty to p rovide additional constraints 
on initial conditions [29-|4a|. 



In [47[ the collision energy dependence of particle 
transverse momentum spectra and elliptic flow coeffi- 
cients were studied using (2-|-l)-d ideal hydrodynam- 
ics with longitudinal boost invariance and a bag-model 
equation of state. In this work, we revisit this prob- 
lem using more realistic (2-|-l)-d viscous hydrodynamics 
coupled with a modern lattice QCD based equation of 
state [13, We also study the differences between 

the two most popular initial conditions obtained from 
Monte Carlo versions of the Glauber (MC-Glauber) and 
Kharzeev-Levin-Nardi (MC-KLN) models 

Our work has some limitations which must be kept in 
mind before comparing them with experimental data. As 
the collision energy decreases, the Bjorken assumption of 
longitudinal boost invariance will gradually break down 
[5^ . Furthermore, since the fireball will spend less time 
in the QGP phase, the hadronic phase becomes more im- 
portant and occupies a larger part in its dynamical his- 
tory. Curing these two major shortcomings will require 
(3+l)-d viscous hydrodynamic simulations 12^ . [5^ cou- 




pled with a microscopic hadronic afterburner 
The present work does not aim at extracting precise in- 
formation of QGP transport properties from a compari- 
son with experimental data. Its main purpose is to ex- 
pose systematic quantitative trends in observables as a 
function of collision energy in the relativistic heavy-ion 
collisions. 

In the next section, we will describe the setup of our 
models and discuss our parametrization of the initial con- 
ditions as a function of -y/s. In Sec HI, we will present 
trends for the transverse momentum spectra and differ- 
ential elliptic flow for charged hadrons as increases 
from 7.7AGeV to 2760 A GeV. Identified particle spec- 
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tra and their elliptic flow V2 will be discussed in Sec. IV. 
In Sec. V, we generalize the definition of the spatial ec- 
centricity to an isothermal hyper-surface. Based on this 
generalized formulation, we perform a shape analysis on 
the final kinetic freeze-out surface and study the depen- 
dence of the final eccentricity on y^. Sec. VI is devoted 
to some concluding remarks. 



II. EVOLUTION OF CHARGED HADRON 
MULTIPLICITY AND TOTAL ELLIPTIC FLOW 

In this work, we employ the (2-|-l)-dimensional viscous 
hydrodynamic model VISH2+1 which im plem ents boost- 
invariance in the l ong itudinal direction |14{ . Similar to 
past work SJ HI, IH, [H], we use two different 

types of initializations taken from the MC-Glauber and 
MC-KLN models to generate initial entropy density pro- 
files. Over one million Monte Carlo events are gener- 
ated and sorted into collision centrality bins according to 
their number of participant nucleons. Each event is re- 
centered to the beam axis and rotated in the transverse 
plane such that its minor axis aligns with the impact pa- 
rameter. Then we average the events to obtain a smooth 
average initial entropy density for each centrality bin. 

In the MC-Glauber runs we take for the specific shear 
viscosity the value 77/s = 0.08 since this value was shown 
in m, m, Hll, 11^ to provide a reasonable description of 
the charged hadron V2(j>t) and vsipr) data measured 
by the RHIC and the LHC experiments For 
the initial entropy density we make the two-component 
ansatz 



1 



(1) 



with a wounded nucleon (WN) to binary collision (BC) 
mixing ratio x that is adjusted to reproduce the mea- 
sured centrality dependence of the final charged hadron 
multiplicity density dNch/drj. For Pb-|-Pb collisions at 

= 2760 A GeV we use a; = 0.118 as determined in [ilj . 
For Au-|-Au collisions at RHIC energies, we determine n 
and X in Eq. (1T|) by a two-parameter fit to the RHIC 
data at ^ 200 AGcV [591-163} obtaining x = 0.14. For 
RHIC collisions at lower y/s we keep the mixing ratio 
fixecfl at a; = 0.14, tuning only the normalization fac- 
tor K to reproduce the charge multiplicity in the 0-5% 
most central collisions. For = 63 A GeV, the desired 
charged multiplicity is taken from experiment [soj . For 

< 63 A GcV, we presently lack experimental informa- 



^ The main reason of keeping the mixing ratio fixed is because the 
measured centrahty dependence of charged muhipHcity for the 
lower energy runs has not yet been pubUshed. We do not expect 
quaUtative changes to the conclusions drawn in this paper once 
our assumed values will be replaced by actual measurements. 



(A GeV) 


To (MeV) 


Tf— To (fm/c) 


dNch/dr) 


AuAuO 7.7 


269.2/233.7 


9.3/9.1 


212.3/212.1 


AuAu@ 11.5 


287.5/252.0 


10.0/9.8 


266.7/266.4 


AuAu@ 17.7 


304.8/269.8 


10.5/10.3 


325.3/324.9 


AuAu@ 19.6 


308.7/274.3 


10.6/10.4 


339.2/338.8 


AuAu@ 27 


320.1/286.4 


10.9/10.7 


382.9/382.1 


AuAu@ 39 


332.2/298.9 


11.2/11.0 


432.7/432.3 


AuAu@ 63 


341.1/306.4 


11.4/11.2 


472.0/472.9 


AuAuQ 200 


378.6/347.0 


12.2/12.1 


661.9/690.0 


PbPb@ 2760 


485.2/443.9 


14.2/14.2 


1575.7/1597.2 



TABLE I. The initial temperature at the center of the fireball, 
fireball lifetime and final charged hadron multiplicity of 0-5% 
most central collisions are listed. The results on the left are 
from MC-Glauber initial conditions with r]/s = 0.08, the right 
are for MC-KLN with ry/s = 0.2. 



tion and therefore use the empirical formula [47 



drj 



312.5 logioA/s- 64 



(2) 



The actually employed final charged multiplicities are 
listed in Table H] 

The MC-KLN calculations are done using a Monte- 
Carlo sample of initial profiles with identical properties 
as those used in [2^. These initial MC-KLN profiles 
were evolved hydrodynamically with a larger viscosity 
77/s = 0.2 to compensate for the larger initial eccentrici- 
ties. Again, this choice was shown to yield a good overall 
description of the measured transverse momentum spec- 
tra and elliptic flow in 200 A GeV An- An collisions at 
RHIC [ill and gave an impressively accurate prediction 
[2^ [131 tti'^ unidentified and identified charged hadron 
spectra and ellip tic fiows in 2760 A GcV Pb-Pb collisions 
at the LHC [il[33l. The large rj/s of 0.2 fails, however, 
to reproduce the large W3 measured in Pb+Pb collisions 
at the LHC [3l[4l|. In the MC-KLN model, the initially 
produced gluon density profile is controlled by the de- 
pendence of the saturation scale Qs on the position x± 
in the transverse plane. For a nucleus with mass number 
A it is given by [52| 



Ta{x±) (Xq 

Ta,o 



(3) 



where Ta{x±) is the nuclear thickness function. We use 
the same parameter set((5^g = 2GeV^, Ta^o = 1.53 fm^, 
A = 0.28 and xq = 0.01) as proposed in Ref. [H]. For 
Au+Au at 200 A GeV and Ph+Ph at 2760 A GeV the nor- 
malization constant for the initial entropy density was 
determined by an overall fit to the centrality dependence 
of These best fits result in slightly different 

values for the 0-5% most central collisions than obtained 
for the corresponding MC-Glauber cases (sec Table|l]). At 
lower energies, the normalization factor was again fixed 
to reproduce the desired charged hadron multiplicity den- 
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sity dNch/dr] for the 0-5% most central collisions for all 
(see Table H. 

Starting from an assumed thermalization time tq = 0.6 
fm/c [23,!2i,[4l|, both the MC-Glauber and MC-KLN ini- 
tial conditions are evolved hydrodynamically using the 
lattice QCD based equation of state (EOS) s95p-PCE 
[27I lisj . This EOS accounts for chemical freeze-out at 
Tchom = 165 MeV before thermal decoupling which is 
taken to occur along an isothermal surface of temper- 
ature Tdoc = 120 MeV. We convert the hydrodynamic 
output along the kinetic decoupling surface into final 
hadron distributions using the Cooper-Frye prescription 
[65| . Strong resonance decays are taken into account up 
to 2 GeV in particle mass. 

We point out that we keep the value of the specific 
shear viscosity 77/s unchanged as we go to lower collision 
energies. As the highly viscous hadronic phase becomes 
more and more important at lower collision energies, vis- 
cous hydrodynamic simulations with temperature inde- 
pendent 77/s will eventually break down. Worse, larger 
77/s values in the hadronic phase jeopardize the validity 
of the viscous hydrodynamic approach altogether (66j . In 
this study, we arc not trying to extract the temperature 
dependence of 77/s from a serious comparison with exper- 
imental data; our goal is to present a systematic study 
of the -y/s dependence of hydrodynamic variables. For 
this reason, we run viscous hydro all the way down to 
y/s ~ 7. 7^ GeV with constant ij/s compared to [l^l, our 
simulations are more realistic by including viscous effects 
in the hydrodynamic evolution and using a better EOSQ 
Also, we here study two different initialization models 
and include (at least on average) the effects of event- 
by-event fluctuations whereas in [43] an optical Glauber 
model was used for initialization which gives too small 
eccentricities in most central collisions. We will see that 
the different ^/s dependence from the two initialization 
models will help us to further distinguish between the 
two initialization models. 

In Tabic |I] we have summarized the global variables 
for our hydrodynamic simulations. At higher collision 
energies the evolution starts with a higher peak initial 
temperature, thus probing the nuclear matter at higher 
temperature and resulting in a longer lifetime of the fire- 



^ It should be noted that our EOS assumes zero net baryon density 
— an assumption that is untenable in the lower half of the colli- 
sion energy range explored here. To include effects from non-zero 
baryon density would require an upgrade of VISH2+1 to solve ad- 
ditionally for the space-time evolution of the conserved baryon 
current. This is important for the correct prediction of the final 
baryon and meson abundances at lower y/s which our present 
code can not achieve. However, what matters for the evolution 
of radial and elliptic flow is the stiffness of the EOS, embodied by 
the pressure (whose gradients supply the hydrodynamic acceler- 
ation) and its relation to the energy density (inertia) of the fluid, 
p(e, n). Since, for not too large baryon densities n, this relation 
depends on n only very weakly [gJI, the use of a baryon- free EOS 
is expected to work well for the systematic flow study presented 
here. 



ball. At LHC energy we find a peak temperature that 
is about twice as large as that reached at the lowest 
collision energies at RHIC, and the lifetime is about 5 
fm/c longer. MC-Glauber initial conditions have about 
30 MeV higher peak temperatures than MC-KLN ones. 
This is mostly due to the fact that the specific shear vis- 
cosity in the MC-KLN runs is about 2.5 times larger than 
for MC-Glauber runs, causing stronger viscous heating 
and larger entropy production during the hydrodynamic 
evolution. The same final multiplicity can thus be 
reached starting from less initial entropy. A larger spe- 
cific shear viscosity also helps the system to develop more 
radial fiow in the transverse plane, by speeding up the 
equalization between transverse and longitudinal veloc- 
ity gradients (the latter are initially very large). This 
larger transverse expansion rate compensates for the vis- 
cous heating effects on the lifetime, resulting in a slightly 
shorter lifetime for the MC-KLN runs. 

In Figs. [lja,c) we show the centrality dependence of 
the charged hadron multiplicity for both MC-Glauber 
and MC-KLN models with collision energies from -y/s = 
7.7 to 2760 A GeV. The reader should note that all results 
in Fig. [1] account for viscous entropy production during 
the hydrodynamic evolution. We checked that at LHC 
and top RHIC energies (top two curves in Figs. l(a,c)) 
our results for both initialization models agree well with 
the experimental data j59l]63l IH, . Our lower collision 
energy predictions can in the future be checked against 
data collected in the RHIC BES program. 

In order to study how the centrality dependence 
changes with -y/s, we scale in Figs. ^p,d the lower col- 
lision energy results by constant factors to align them 
with the LHC curve in central (0-10%) collisions. For the 
MC-Glauber model we find good -yi-scaling: the curves 
almost fall on top of each other. For the low energy runs 
at RHIC this is, of course, sensitive to the fact that we 
keep the mixing ratio between the wounded nucleons and 
binary collisions fixed, and it also reflects the fact that 
viscous entropy production is small and has little effect 
on the centrality dependence. On the other hand, for 
the MC-KLN model the slope of the centrality depen- 
dence gets flatter as the collision energy decreases. Only 
the top RHIC and LHC energy curves approximately fall 
on top of each other; at lower energy this y^-scaling is 
broken. We found that this tendency originates in the 
nature of the MC-KLN model itself: Even though vis- 
cous entropy production is larger (due to the larger 77/s 
used in the MC-KLN runs) , its centrality dependence has 
only a minor effect on the centrality dependence of 
and cannot explain the different shapes of the curves in 
Figs. [T]D,d. Our MC-KLN calculations thus predict a vi- 
olation of the -^/s-scaling of the centrality dependence of 
at lower collision energies that is not seen with the 
MC-Glauber initial conditions. This may help to dis- 
criminate experimentally between these models. 

A "universal" scaling behavior of the eccentricity- 
scaled elliptic flow as a function of charged hadron mul- 
tiplicity density ("multiplicity scaling") [t^ was stud- 
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FIG. 1. (a): Centra- 
lity dependence of final 
charged hadron multipli- 
city per participant nu- 
cleon pair as a function 
of TVpart for MC-Glauber 
initial conditions, with 
collision energies vary- 
ing from -^s = 7.7j4 GeV 
to ^ = 2760AGeV. (b): 
Centrality dependence of 
i^sh from the lower en- 

(177 

ergy runs in (a) scaled up 
to the LHC results, for 
shape comparison. (c, 
d): Same as (a, b) but 
for MC-KLN initial con- 
ditions. 
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FIG. 2. Eccentricity- 
scaled pT-integrated U2 
plotted as a function of 
the charged hadron mul- 
tiplicity density for differ- 
ent collision energies, for 
MC-Glauber initial con- 
ditions with Tj/s = 0.08 
(a) and MC-KLN profiles 
with 77/s = 0.2 (b), re- 
spectively. 



ied within viscous hydrodynamics in jlSj and was later 
used to extract the specific shear viscosity from -y/s = 
200 A GeV Au+Au collisions at RHIC [Hi. The authors 
of [2^ and [m found that this "universal" scaling breaks 
down as y/s increases but disagreed on the sign of the 
sealing breaking effects. In Fig. [2] we explore the break- 
ing of "multiplicity scaling" over a wider range of -^/s, for 
both of the initialization models. For MC-Glauber ini- 
tial conditions (Fig. [2ji) eccentricity-scaled elliptic flow 
shows surprisingly good universality of the "multiplic- 
ity scaling" curve as the collision energy varies from 7.7 
to 2760 j4 GeV: The curves at different -^s fall almost 
perfectly on top of each other. For MC-KLN (Fig. Eb), 
on the other hand, the "universal scaling" breaks in the 



same direction as previously shown in |28j : lower colli- 
sion energies result in larger V2/t2 values at the same 
charged hadron multiplicity density. We found that the 
main reason for the different collision energy dependence 
between the MC-Glauber and MC-KLN models lies in 
the different centrality dependences of the initial overlap 
area in the two models. The initial overlap area is cal- 
culated as 5 = 7r^(a;2)(2/2)^ where (x^) = J^i^^.^].'^ 
is evaluated with the initial energy density as weight 
functionlfl As the collisions become more peripheral, 



^ The initial entropy density can also be used as weight. In [28ll 
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FIG. 3. (a,c): Transverse 
momentum spectra of all 
charged hadrons from 
central (0-5% centrality) 
Au+Au and Pb+Pb col- 
lisions at 0-5% centrality 
for different collision 
energies. (b,d): the cor- 
responding differential 
elhptic flow at 20-30% 
centrality. 
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the overlap area in the MC-KLN model decreases more 
rapidly than in the MC-Glauber model. In Pb-|-Pb col- 
lisions at = 2760 A GeV, the overlap area S for 
MC-KLN decreases from 23.6 fm^ in the 0-5% most cen- 
tral collisions to 4.7 fm^ in the 60-70% centrality class; 
for MC-Glauber, S decreases instead from 22.8 fm^ to 
6.5 fm^. This slightly faster drop of the overlap area in 
the MC-KLN model shifts the "universal" scaling curves 
in Fig. [2] to the right and shrinks the covered range in 
{\/ S)dNch/ drj. We further checked that the centrality 
dependence of the overlap area changes little as y^s varies 
from 7.7 to 2760^ GeV. The different ys-dcpendences 
of V2/e2 as function of dNch/drj in Figs, [^i andl^b thus 
reflect primarily the fact that the shape of the initial 
profiles evolves differently with centrality in the two ini- 
tialization models. Fig. [5] can thus be used to check ex- 
perimentally the consistency of the centrality dependence 
of the source size and shape in the initialization models. 

III. CHARGED PARTICLE pt-SPECTRA AND 
DIFFERENTIAL ELLIPTIC FLOW 

Figures[3^,c show the pT-spectra of all charged hadrons 
in the 0-5% most central collisions. For both the MC- 



wc showed that the scaling breaking behavior is independent of 
the choice of weight function. 



Glauber and MC-KLN models the slopes of the pT- 
spcctra get flatter as ^/s increases: At higher collision 
energy the fireball lifetime is longer, which allows the sys- 
tem to develop more radial flow. The additional radial 
flow pushes more particles into the high-p^ region, thus 
flattening the spectra. From ^/s = 7.7 to 2760 A GeV, 
the mean pt, (pt> = I dprPTj^/ J dpT^^^g^, in- 
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FIG. 4. Evolution with ^/s of the average transverse mo- 
mentum of charged hadrons from central (0-5% centrality) 
Au+Au and Pb-|-Pb collisions, for MC-Glauber and MC-KLN 
initial conditions. 
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FIG. 5. Evolution with 
-^s of the difTerential 
charged hadron elliptic 
flow V2^{pt,\/s) at 5 
flxed pt values. 



creases by 43% (from 0.48 to 0.68GeV/c) for the MC- 
KLN model and by 140% (from 0.29 to 0.63GeV/c) for 
the MC-Glauber model (see Fig.g]). 

The differential charged hadron elliptic flow is shown 
in Figsl3b,d, for 20-30% centrality. With MC-Glauber 
initial conditions the differential elliptic flow for pT < 2 
GeV remains almost unchanged for y^s > 39 ^4 GeV. Be- 
low 39 A GeV the slope of V2 (pr) begins to decrease. This 
tendency is indeed observed in the RHIC BES experi- 
ments [3, 0, We emphasize that our EOS s95p-PCE 
has no pronounced soft point in the phase transition re- 
gion. This means that the often highlighted "satura- 
tion" of V2{pt) above ^/s = 39A GeV can not be as- 
sociated with a softest point in the transition region. 
It is rather caused by a subtle cancellation of opposite 
yi-dependences of the differential Vs{pt) from light and 
heavy particles (see Fig. |8] below). 

For the MC-KLN model, the slope of the differen- 
tial V2{pt) decreases monotonically and continuously 
with decreasing collision energy. For a temperature- 
independent specific shear viscosity, r]/s — 0.2, the col- 
lision energy dependence of the differential elliptic flow 
observed here is somewhat inconsistent with the exper- 
imental observation of a v'^{pt) that docs not change 
between = 39 and 2760^ GeV. Within the MC-KLN 
framework, this might be taken as an indication for a pos- 
sible temperature dependence of 77/5 [2^[6l,[7l|. Addi- 
tional studies are, however, necessary to fully address this 
issue [Hi. For = 7.7 and 11.5^ GeV, the differential 
V2 is seen to increase more quickly above pr > 2.5 GeV. 
We find this to be caused by large 5f corrections (i.e. 
non-equilibrium corrections arising from non-zero shear 
stresses at freeze out [l^). The larger 5f corrections at 
lower collision energies indicate a narrowing of the tem- 
poral interval during which viscous hydrodynamics is a 
valid description. At lower pT [pt < 2 GeV) , our results 
show monotonic ^/s dependence. 

To further illustrate this point we plot in Fig. [5] the 
•yi-dependence of V2^{pt) at 5 fixed px points. In this 
representation one sees that for the MC-Glauber model 
at any fixed p^-value features as a function of ^/s a 
very broad maximum somewhere around top RHIC en- 
ergy (200 A GeV); for low pr < 0.5GeV/c, this maxi- 



mum occurs at lower ^/s. A similar behavior was seen 
in for ideal hydrodynamics with a bag-model equa- 
tion of state which features a strong minimum ( "softest 
point") in the speed of sound at the phase transition 
temperature. We see that the existence of this maximum 
does not depend on the appearance of a softest point 
in the EOS. Compared to the earlier ideal fluid calcu- 
lations, the position where at fixed pT assumes its 
largest value has been shifted to larger ^/s values by vis- 
cous effects. This shift is seen to be even stronger in the 
MC-KLN case (Fig. ^jp) where the fluid is much more 
viscous. For shear viscosities as large as those needed 
to describe the vf^ measured in 200 A GeV Au-|-Au col- 
lisions with MC-KLN initial conditions (77/s ~ 2.5/(47r)) 
[13, [3 J '^2^{pt,\^) at fixed pr has not yet reached its 
maximum value even at top LHC energies (except for 
very small px < 200MeV/c). 

At the lower end of the -y/s-rangc studied in Fig. O 
the increase with collision energy of V2{pt, \fs) at fixed 
Pt is a consequence of increasing fireball lifetimes which 
allow the initial spatial eccentricity of the fireball to con- 
vert more fully into anisotropic hydrodynamic fiow. At 
higher collision energies eventually the point is reached 
where this momentum anisotropy is fully saturated be- 
fore the system falls apart; longer fireball lifetimes will 
then no longer lead to more anisotropic flow, only to 
more radial fiow. Stronger radial flow, however, pushes 
the momentum anisotropy out to larger pT, by generat- 
ing flatter px distributions. As a result, elliptic flow at 
fixed Pt begins to decrease. In practice, this radial flow 
driven decrease of V2{pt) at fixed px sets in even before 
the px-intcgratcd total charged hadron elliptic fiow vf^ 
has reached saturation [t^, and it accelerates thereafter. 



IV. pT-SPECTRA AND ELLIPTIC FLOW OF 
IDENTIFIED HADRONS 

We now proceed to study how hydrodynamical flow 
affects identified particles. 

It is well known that thermal spectra from a 
static fireball exhibit mT-scaling, dNi/ (2ndymTdmT) ^ 
y/rriT e~™^/-^ [t^, and that radial fiow breaks this seal- 
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ing. In order to isolate the radial flow effects we therefore 
plot in Fig. [S] the mr-spectra of identified particles as a 
function of mT—rriQ, for four selected \/s values. Ex- 
cept for minor effects from the viscous 5f corrections, 
resonance feed-down and Bose statistics for pions, in 
the absence of flow the slopes of the rriT-spcctra would 
be the same for all liadron species. To show the flow- 
induced slope difference, we scaled in Figs. |6] the heavy 
particle spectra by constant factors to the same value at 
TOy— "T'O = 3 GeV. At large mT—mo rest mass effects be- 
come negligible, and all hadrons have approximately the 
same inverse slope TeS = Tdocy/{l+{v±))/{l-{v±)) [zl]. 
In Fig. [S] we find that for low -^5 values rriT-scaling is sig- 
nificantly broken only at rriT—mo < 2 GeV, while at LHC 
energy the flow-induced breaking of my-scaling extends 
to 3 GeV of transverse kinetic energy. At low toj^— toq 
the spectra are split by hadron mass effects, and this 
splitting increases with ^/s due to the increasing radial 
flow which pushes heavier particles to larger pT- At the 
highest collision energy ^/s = 2760 A GeV we observe a 
particularly strong concavity of the pion spectra at low 
niT—mQ^ due to Bose statistics. For other mesons, their 
heavy rest masses suppress Bose effects. 

The flow-induced breaking of mr-scaling is seen even 
more clearly when one plots heavy-to-light particle ra- 
tios (such as p/tt^, A/K~^) as a function of transverse 
kinetic energy. For a static thermalizcd fireball, these ra- 



tios should be independent of mr— mo, up to small quan- 
tum statistical corrections arising from the pion spectra 
at small m^— mp. Fig. [7] shows that for an expanding 
fireball these ratios increase with increasing transverse 
kinetic energy, at a rate that itself increases with -y/s, 
reflecting the larger radial flow at higher collision ener- 
gies. A little more careful inspection and thought reveal 
that, in fact, stronger radial flow increases the p/tt"*" and 
A/ ratios at large mT—mQ while decreasing them at 
small mT—rriQ. This is so because in our simulations the 
PT-integrated particle ratios are the same at all collision 
energies, as we assumed zero baryon chemical potential 
and the same chemical and kinetic freeze-out temper- 
atures at all y/s. In addition, radial flow flattens the 
TOT-dependence of these ratios at low rriT—rno, due to 
the "flow shoulder" developing in the heavy-particle m^- 
spectra at low transverse kinetic energy when radial flow 
gets strong. This shoulder is weaker for protons than for 
A's, but in the p/ir^ ratio the more prominent Bose effect 
in the pion spectra at high collision energies additionally 
helps to flatten out the p/n'^ ratio at small mT—mQ. 
Overall, Fig. [7] shows that these features are all very sim- 
ilar for MC-Glauber and MC-KLN initial conditions. 

In Fig. [8] we show the differential elliptic flow of 7r+ 
and p for ^/s = 7.7 to 2760 A GeV. For pions, the differen- 
tial V2 varies with yG very similarly to the total charged 
hadron elliptic flow shown in Figs. [3]3,d. For protons the 
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and A/K'^ (b,d) ratios as 
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Please note the dramatic 
increase of radial flow 
effects on these ratios 
between RHIC and LHC 
energies. 



strong radial flow "blueshifts" the entire elliptic flow to 
higher pT- So for higher collision energies, the values of V2 
are smaller in the low p^-region and larger in the high-py 
region. We flnd that below 200 A GeV the proton V2{pt) 
at low pt is almost independent of ^/s. At LHC en- 
ergy, on the other hand, the blucshift is really dramatic, 
reflecting the much stronger radial flow at this high colli- 
sion energy. The total charged hadron elliptic flow is the 
combination of contributions from light pions and less 
abundant heavy particles. Since with increasing collision 
energy the elliptic flow of heavy particles decreases at 
low PT, they effectively cancel the weak increase of the 
light pion V2- This results in the apparent saturation of 
charged hadron differential elliptic flow over a wide pt 
range from ^/s = 39-2760 A GcV that was seen in Fig.jSjj 
for MC-Glauber model. For the MC-KLN model, this 
cancellation is less efficient because the increase with -^/s 
of the pion V2 is stronger (see Fig. |8l;,d). Therefore, for 
the MC-KLN model the total charged hadron elliptic flow 
keeps increasing as the collision energy increases. 



time T by 



£x{t) 



J dx dy (y^-x^)7e(x,^;T) 
/ dx dy (y2+a;2)7e(x, y; r) ' 



(4) 



The weight function 7e(x, y; t) is the energy density in 
the laboratory frameo 

Since the measured hadrons are only emitted from the 
final kinetic freeze-out surface, experimentalists can only 
infer the shape of that surface, by exploiting two-particle 
momentum correlations among the emitted particles and 
their dependence on the azimuthal angle around the 
beam axis [t^ . [TSj . For comparison with such experi- 
mentally determined final source eccentricities [76l479| . a 
more meaningful theoretical quantity would be the spa- 
tial eccentricity of the final freeze-out surface E, 



ut^d^afj_ (y'^+x^) ' 



(5) 



instead of a constant proper time surface. In Eq. ^ 
is the flow velocity on the surface S, and we used the fact 



V. SPATIAL ECCENTRICITY AT FREEZE-OUT 

We conclude this paper by presenting a novel shape 
analysis of the evolving flreball. Theoretically, the spatial 
eccentricity is conventionally dcflncd at flxcd proper 



* The entropy density s{x, y; t) can also be used as weig ht function 
for calculating the eccentricity. The authors of [Sal concluded 
from a study of fluctuating initial conditions that these two defi- 
nitions yield initial spatial eccentricities that arc linearly related 
to each other although the actual values are slightly different. 
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that for our EOS an isothermal freeze-out surface is also 
a surface of constant local energy or entropy density, and 
that therefore the weight functions e(x, y, r) or s(a;, y, r) 
cancel between numerator and denominator^ Thus the 
spatial eccentricity ([5]) defined on an isothermal surface 
is independent of whether we weight with energy or 
entropy density. 

In Fig. [51 we show the final eccentricity calculated 
along the kinetic freeze-out surface, Tdoc = 120 MeV, 
as a function of collision energy. For both MC-Glauber 
and MC-KLN models, as the collision energy increases, 
the final spatial eccentricity ef decreases monotonically. 
This is because at higher collision energy the system lives 
longer, giving the fireball more time to decompress and 
(due to anisotropic fiow) become less deformed. For suf- 
ficiently large initial energy density, the fireball has actu- 
ally enough time to become elongated along the reaction 
plane, instead of its original elongation perpendicular to 
it [s^iUll- In Fig- Hlwe compare our results with recent 
(preliminary) STAR data from an azimuthal HBT analy- 
sis [Sli. MC-Glaubcr runs with 77/5 = 0.68 quantitatively 
reproduce the data at ^fs = 200^ GeV while underprc- 



^ If we do not cancel the weight function and replace E by a con- 
stant proper time surface, the definition JSj reduces to Eq. l|4]l 
for longitudinally boost-invariant systems. 



dieting the final eccentricity by ~10% at lower energies. 
MC-KLN initial conditions with 77/5 = 6.2 result in 15- 
20% larger final eccentricities than both the MC-Glauber 
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FIG. 9. v^-dependence of the final spatial eccentricity ei of 
the isothermal kinetic freeze-out surface at Tdoc = 120 MeV, 
for 10-30% centrality. The initial eccentricity is 0.26 for the 
MC-Glauber model and 0.32 for the MC-KLN model. The 
experimental points indicate preliminary data [8l[ from an 
azimuthal HBT analysis by the STAR Collaboration. 
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runs and the STAR data, due to the ^20% larger initial 
eccentricities of the MC-KLN profiles. In Table |T] we see 
that, at the same ^/s, the fireball lifetimes with MC-KLN 
and MC-Glauber initial conditions are very similar. In 
spite of the faster evolution of radial flow for the MC- 
KLN initial conditions, the larger initial eccentricity in 
the MC-KLN model is preserved all the way to the end 
of the hydrodynamic evolution. Extending our calcula- 
tions to LHC energy we predict that et will approach 
zero around ^/s = 2.76-5.5 A TeV. Again, this is a result 
of the longer fireball lifetime at LHC energies. For even 
larger -^s, ef will turn negative, i n q ualitative agreement 
with previous calculations in |82l . |83| using ideal fluid dy- 
namics and a less realistic EOS. Contrary to our present 
work, the authors of [lil, H^l calculated the azimuthal 
HBT radii from the Cooper-Frye output of their hydro- 
dynamic simulations. In future work we will calibrate our 
definition ([5]) for Sf against the final eccentricity value ex- 
tracted from azimuthal oscillations of HBT radii. Here 
we only note that an azimuthal HBT analysis at LHC en- 
ergy will help to further test predictions from the viscous 
hydrodynamic model. 

The careful reader may have noticed that there is a 
small kink at ^/s = 63AGeV in the slopes of the solid 
and dashed lines in Fig. [9] and earlier in Figs. |4] and 
[SJ Unfortunately, this is not a phase transition signa- 
ture but rather an artifact from our normalization of the 
initial energy density profile to experimentally measured 
final charged multiplicity data for y/s > 63 A GeV and 
to the empirical formula Eq. ^ for < 39 A GeV. At 
-y/s = 63 A GeV the experimentally measured dNch/drj is 
~ 5% smaller than the value obtained from Eq. ([2]). If we 
use Eq. ([2]) instead of the measured value to normalize 
our initial energy density profiles at i/i = 63A GeV, the 
fireball lifetime increases slightly, decreasing the final ec- 
centricity by a few precent and removing the kink in the 
theoretical curves. The origin of an apparently similar 
kink in the STAR data which happens to occur (we be- 
lieve: accidentally) also around ^/s = 63^ GeV deserves 
further study. 



VI. SUMMARY 

Using (2+l)-dimensional viscous hydrodynamics with 
properly implemented hadronic chemical freeze-out at 
Tchcm = 165 MeV, we have studied systematically the evo- 
lution of hydrodynamic observables with collision energy 
in the range 7.7 < -/s < 2760 ^ GeV. Over this range of 
energies, the initial peak temperature almost doubles and 
the fireball lifetime increases by about 60%. We find that 
for temperature independent specific shear viscosity the 
MC-Glauber model shows almost perfect "multiplicity- 
scaling" of the eccentricity-scaled charged hadron ellip- 
tic flow. For the MC-KLN model this scaling is broken: 
as y/s increases, the v'^je-i vs. {l/S){dNch/dy) curves 
shift to the right (Fig. ^jp). We found that this breaking 
of multiplicity scaling in the MC-KLN model originates 



from a steeper centrality dependence of the nuclear over- 
lap area. 

For both initialization models, higher collision ener- 
gies generate stronger radial flow which results in flat- 
ter hadron spectra and a corresponding increase of the 
mean px of charged hadrons. For the MC-Glaubcr model 
we observed an approximate "saturation" of the charged 
hadron differential elliptic flow at fixed px in the region 
\/s > 39 ^ GeV, similar to what is observed experimen- 
tally. We believe, however, that the word "saturation" 
describes the observations incorrectly and that what is 
seen is better described as a very broad maximum (as 
a function of y/s) of the differential elliptic fiow at fixed 
Pt, caused by the interplay of (i) growing total momen- 
tum anisotropy (which increases V2) and increasing ra- 
dial flow (which decreases V2 at fixed px by shifting it 
to larger px), and (ii) increasing V2{pt) for pions and 
decreasing V2{pt) for kaons, protons and other heavy 
hadrons between RHIC and LHC energies. The mech- 
anism (i) causes maxima of V2{pt) at fixed pt for all 
hadron species, but located at lower -y/s values for heav- 
ier than for lighter hadrons (due to the mass-dependence 
of radial flow effects on the pT-spectra) . The mechanism 
(ii) ensures that the maximum for all charged hadrons 
is broader in y/s than for each hadron species individu- 
ally and thus manifests itself as a broad plateau that (for 
7]/s = 0.08) happens to span the collision energy range 
from upper RHIC to LHC energies. The position in ^/s 
of the maximum of V2(j>t) at fixed pt for each hadron 
species depends on the viscosity of the fluid (which con- 
trols the interplay in the development of radial and ellip- 
tic flow during the fireball expansion) and increases with 
77/s. For MC-KLN initial conditions, which require ^2.5 
times larger 77/s for a successful description of elliptic 
flow data at RHIC and LHC, v'^[pt) at fixed px has not 
yet reached its maximal value even at LHC energies. 

Finally, we have proposed an improved measure for the 
final fireball eccentricity at kinetic freeze-out and studied 
its evolution with collision energy. It is found to decrease 
monotonically with increasing collision energy, at a rate 
that is roughly consistent with recent experimental mea- 
surements. Its absolute value agrees with the data better 
for the MC-Glauber than for the MC-KLN model - the 
~20% larger initial eccentricities of the MC-KLN profiles 
yield final freeze-out eccentricities that again appear to 
be ^--^20% larger than those from MC-Glauber initial pro- 
files, and lie significantly above the measured values. Nei- 
ther model describes the available data perfectly; in view 
of the limitations of the purely hydrodynamic approach 
employed here (c/. our discussion in the Introduction) 
this is not too surprising. The model predicts, however, 
robustly that at top RHIC energies the final freeze-out 
source is still out-of-plane elongated (as experimentally 
observed), but that at LHC energies the final eccentric- 
ity should approach zero. Measurements that test this 
prediction should soon become available. 
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